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I.  INTRODUCTION 

Recently  a technique  has  been  developed  [1]  for  extending  the  use 
of  the  Geometrical  Theory  of  Diffraction  (GTD)  by  using  the  method  of 
moments  to  represent  the  current  near  a discontinuity  in  a surface  such 
as  near  the  edge  of  a two  dimensional  wedge.  The  wedge,  of  course,  is 
a canonical  problem  in  the  GTD  and  is  used  merely  as  a basic  geometry 
to  test  the  method. 

The  previous  work  [1]  in  extending  the  GTD  via  the  moment  method 
required  special  consideration  for  those  cases  where  the  incident  ray 
was  nearly  tangent  to  one  of  the  faces  of  the  wedge  and  the  results  were 
somewhat  dependent  upon  the  location  of  the  match  points  used.  In  this 
report  these  difficulties  are  overcome  by  using  a series  of  three  terms 
based  on  approximate  expressions  for  the  Fresnel  integral. 

Results  for  the  currents  on  the  faces  of  the  wedge  will  be  shown 
for  a variety  of  cases.  In  all  cases  the  currents  obtained  by  the 
improved  combined  GTD  moment  method  formulation  are  compared  with  the 
exact  currents  and  the  agreement  is  found  to  be  excellent. 


II.  FORMULATION 

Consider  the  general  wedge  diffraction  situation  shown  in  Figure 
1.  The  total  magnetic  field  at  |ny  point  of  the  space  is  gjven  as  the 
summation  of  the  incident  field  H’  and  the  scattered  field  Hs.  At  the 
point  P(p-|  ,$i  ,0i ) the  field  will  be 


H = H1  + 


(1) 


The  scattered  field  is  expressed  in  terms  of  the  vector  potential  A by 

ti5  - 1 ? x a (2) 

where  ^ is  related  to  the  wedge  current  by  the  following  expression 


^(p-j)  = n JJ  3(r)  g(p1 ,r)ds' 


(3) 


sw 


(sw  is  the  wedge  surface  and  G is  the  Green's  function). 


If  we  suppose  that  the  wedge  has  perfect  conductivity  then 


H 


Figure  1.  Wedge  diffraction  geometry  used  in  MM-GTD  solution. 


J(r)  = nr  x H(r) 


where  nr  is  a unit  vector  normal  to  the  wall.  By  the  use  of  the  above 
equations  we  can  find  that: 


J(r)  = n x M1  - n x 3(r')  x 7 G(r,r’)ds’  . 

t J ! ; 


Given  that  the  magnetic  field  as  well  as  the  vector  potential  are 

-jkz  cos 60 

dependent  on  z only  through  the  phase  factor  e .all  problems 

we  will  consider  ar^e  two-dimensional.  The  current  may  be  represented 
by  the  components  Jt  on  xy  plane  and  Jz  on  xz  or  yz  plane  (Figure  2) 
and  will  take  the  form: 


SBE32SS52SZ 


Figure  2.  Current  of  a point  p1  used  in 
magnetic  integral  equation. 


-t  -*  -t  ■+  -jkz  cos 

3(r‘)  = [?z(x,y)  + JT(x,y)]e 


(6) 


By  the  help  of  Equations  (5)  and  (6)  we  have  the  integral  expression 
for  the  Jj  and  Jz  currents  as  follows: 


JT(p)  = np  x hP(p)  -np  x j 


^t(p 1 ) x v G(p ,p 1 ) * d A * 


Jz(p)  = np  x^(P)-npx  j ^z(p’)  x 7p  G(p,p')df  - 

c 

-jkz  np  x z x 3t(p')  • G(p.p')  • d*1 


(7) 


(8) 


The  Green's  function  in  the  two  dimensional  case  reduces  to  the  Hankel 
function  [2]  and  is 


G(p,p')  = Hq(  kT  |p-p ' |) 


(9) 
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In  using  the  GTD-MM  solution  [1]  we  must  first  recognize  the  difference 
in  the  diffraction  near  and  far  from  the  edge  of  the  wedge.  Thus  if 
the  plane  of  incidence  is  in  region  (I),  then  we  have  two  reflected 
rays  (Figure  3a)  Ro»R]«  And  we  can  say  that  the  problem  may  be 
separated  into  two  problems.  In  both  of  them  we  have  to  study  the 
diffraction  from  3 different  rays  as  shown  in  Fiqures  3b, c The  in- 
cident angles  for  these  are: 

xz  plane 


1.  * 

2.  360-4, 

3.  360-4>-2y 


The  diffracted  current  on  the  x axis  will  be  the  summation  of: 


yz  plane 

1 . 360-4>-y 

2.  4>- y 

3.  4>+y 


DX(«J>)  • e 


-jkTrcos4> 


/ k'r  r ('l +'co's4~) 


(10) 


y jkTr  cos 4> 

J2x  = Dx(360-$)e  1 


f • 2 

J e"JT  dr 
l/krr(l+cos4>) 


01) 


Dx(360-*-2y)  e 


jk-j-r  cos  (4>+2y) 


i/kr,r(l+cos(4'+2Y)) 


(12) 


where  the  diffraction  coefficients  Dx  are  assumed  to  be  unknown  and 
dependent  on  each  other  and  the  corresponding  CK  on  the  y-axis. 

A Fresnel  integral  can  be  approximated  as  [3]: 


-»J 


i 

i 


0 


(13) 


e"JT  dT 


*2  r~  » 

yy. 


- x > 1 

TT 


where  aQ,  a^  and  a2  are  known  coefficients.  If  the  three  rays  are 
far  from  the  x axis  then 


Jx  = Jlx  + J2x  + J3x  = DxT 

► x 


If  any  of  the  three  rays  is  near  the  x-axis,  then  either 


7 kyr(  1 +cos4> ) > 1 


~ k-j-r  [1  + cos(<fr+2y)]  > 1 


is  true  when  r is  large.  In  our  study  we  can  say  that  if  either  (16) 
or  (17)  is  not  true,  then  the  distance  r is  0 < r < 3A,  and  thus  that 


where 


r > minfr*,^) 


1 4sineQ(l+c 


rx  = * 

2 4 s i7^30L  1 +cos  (♦+2y  ] 


6 


The  approximation  of  the  total  current  on  x-axis  will  be 


1 


0<x<min(r^,r2) 


Jx+Jx  + ^ Dx  n min(ri»r2)  <x<~ 

{ x x n=-l  x (/x)n  1 1 

where  P(x-x  ) is  a pulse  function  with  ax  weiqht. 

m m 3 

The  same  approximation  can  be  used  for  the  current  on  y-axis 


[ ay  P(y-y  ) 


0<y<min(r>',r^) 


1 "JkTy 

JV/+  l Dn  _ 

y y n=-l  y ^ 


min(^,ry)<y< 


If  the  incident  ray  is  in  region  (II)  (Figure  3a),  then  in  expressions 
(19),  (20)  we  will  not  have  J1  or  Jr  in  the  face  which  is  in  the  shadow 
region  of  the  incident  plane.  By  the  help  of  the  equation  (1)  and  with 
the  same  work  as  Burnside  et  al.  [1]  we  can  find  the  unknown  coefficients 

am*  aX*  °X’  Thus* 

Point  on  x-axis 
1.  PM  - region  (x=xp) 


P + J I + ) D*  I*  = -H1 (x  ) - k 
r L m mn  uy  An  z'  n;  n 


2.  GTD  - region  (x=xn) 


U,  «• ! *>  -• 
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the  expression  for  points  on  the  y-axis  can  be  found  by  changing  the  x 
to  y and  the  y to  x.  In  a general  case  for  finding  the  unknowns 
we  need  k+N+6  equations: 

Taki ng 


(27) 


r = [x2  + y2  - 2xny  cos  y]1/2 


Substituting  the  current  Jy  in  Equation  (8)  one  can  find 
way  the  Jz  current.  A more  applicable  expression  for  J2 
be  found  with  the  help  of  Maxwell's  Equations  from  which 
show  that 


■ -j  77  "P  x 7P  n2(p) 

kT 

By  Equation  (28)  we  have: 

,2,..  x 1 


N 


■W  " HT(xn>  - JuST  l ' 

o n=  I 

• ^ i,  * 


where 


nm 


= 1 


COS  B , ? 

+ T-2  [s,nx  Ll  ' K1  0 


? y -X_yCOSY 

H^kjr)  ^ dy 


6y 


m 


‘ran  ■ V1nY  J [H0(kTr)  - 


H7(kxr)  x -y  cosy 
-XT'  1 • -Hr d> 


fiy, 


m 


<• 


oo 

! 


min(r^,r^) 


-jkTy 

H?(k Tr)  x- 
1 T (<y)1 


y -xny  cosy 


dy 


in  the  same 
current  can 
we  can  easily 


(28) 


(29) 


(30) 


(31) 


(32) 
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I = 


V1nY 


I 

min(r^  .rp 


CHo(kTr) 


H?(kTr) 

'T^F~ 


e xn+y  C0SY 


tf)1 


dy 

(33) 


min(r^,rp 


2.  . y -v  C0SY 

H^(kTr)  dy 


(34) 


I"  I 


min^.rp 


C Ho(kTr) 


H^(kjr)  xp+y  cosy 
~T7r  ^ r 


dy 


(35) 


By  the  help  of  Equations  (8)  or  (28)  we  can  find  the  current 
parallel  to  the  edge.  Both  methods  must  show  that: 


= tan6 


(36) 


which  is  also  true  from  the  principles  of  GTD  for  arbitrary  incidence 
angles  if  the  incident  field  is  such  that  the  problem  is  not  purely 
TE  or  purely  TM. 


III.  EXAMPLES 

Numerical  results  have  been  obtained  for  a variety  of  cases. 

The  first  set  of  results  is  shown  in  Figure  4 for  the  situation  where 
a plane  wave  is  incident  normally  ( 0o=tt/2 ) on  the  edge  of  a 90°  wedge. 
The  agreement  between  the  GTD-MM  results  obtained  using  Equations  (21) 
and  (22)  and  the  exact  solution  [4]  is  seen  to  be  excellent. 

Figures  5 and  6 show  results  for  two  cases  where  a plane  wave  is 
incident  obliquely  on  the  wedge  edge.  Again,  agreement  with  the  exact 
solution  is  excellent  for  both  Jj(x)  and  J2(x).  In  Figures  5 and  6 
as  well  as  4,  it  was  only  necessary  to  use  the  term  in  Equations  (19) 
and  (20)  although  the  same  results  are  obtained  if  all  three  terms  are 
used.  This  is  what  one  would  expect  and  is  in  agreement  with  previous 
work  [1],  However,  to  obtain  the  correct  results  in  Figure  7,  it  was 
necessary  to  use  all  three  terms  since  the  incident  rays  are  near 
grazing  on  one  or  the  other  faces  of  the  wedge. 
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Figure  5a.  Current  perpendicular  to  the  edge  along  x-wall  of 


bb.  Current  parallel  to  the  edge  along  x-wall 
wedge  for  oblique  incidence  case. 


Figure  7a.  Magnitude  of  the  currents  along  x-wall  of  30°  wedge  for 
two  normal  incidence  cases  near  the  shadow  boundaries. 


Figure  7b.  Phase  of  the  current  along  x-wall  of  30°  wedge  for  two 
normal  incidence  cases  near  the  shadow  boundaries. 


T—., 


Figure  8 shows  results  of  a grazing  incidence  case  which  shows 
excellent  agreement  with  the  exact  current  along  the  x wall.  As  in 
Figure  7 retention  of  all  three  terms  will  yield  the  correct  results. 
However,  since  the  incident  ray  is  not  close  to  being  parallel  to  both 
walls  as  in  Figure  7,  it  was  found  that  the  d(_1)  and  D°  terms  could  be 

X X 

set  to  zero  in  Equation  (19)  without  effecting  the  results.  This  is 
indicated  by  Table  1 which  shows  the  minimum  number  of  terms  needed 
for  each  case  considered  in  this  report. 


TABLE 

1 

TERMS  REQUIRED  FOR  VALID  RESULTS.  IN 

ALL  CASES, 

ALL  SIX  TERMS 

CAN  BE 

USED  TO  OBTAIN  THE 

CORRECT 

SOLUTION 

Figure 

D(-l ) 

D<°> 

Dn) 

D(-l) 

D<°> 

D0) 

X 

X 

X 

y 

y 

y 

4 

no 

no 

yes 

no 

no 

yes 

5 

no 

no 

yes 

no 

no 

yes 

6 

no 

no 

yes 

no 

no 

yes 

7 

yes 

yes 

yes 

yes 

yes 

yes 

8 

no 

no 

yes 

yes 

yes 

yes 

9 

yes 

yes 

yes 

no 

no 

yes 

10a 

no 

yes 

yes 

yes 

yes 

yes 

10b 

no 

no 

yes 

no 

no 

yes 

Figure  9 shows  the  current  along  the  x wall  when  the  incident 
ray  is  close  to  grazing  along  that  wall.  In  contrast  to  Figure  8, 
here  the  current  phase  shows  rapid  variation  and  again  agreement  with 
the  exact  solution  is  excellent.  Table  1 shows  that  Dj"')  and 


need  not  be  used  but,  as  in  Figure  8,  if  all  terms  in  the  series  are 
used  the  results  are  unchanged. 


Finally,  Figure  10  shows  the  results  obtained  if  other  combi- 
nations of  the  Dx  and  Dy  coefficients  are  used.  It  is  apparent  that 
if  one  or  more  of  the  four  coefficients  required  for  good  results  in 
Figure  9 is  excluded,  then  incorrect  results  are  obtained.  Moving  the 
match  points  to  different  locations  in  the  GTD  region  does  not  affect 
the  incorrect  results  in  Figure  10  or  the  correct  results  in  the 
previous  cases. 
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Figure  8.  Current  along  x-wall  of  90°  wedge  for  normal  incidence  case 
near  the  shadow  boundary  of  y-wall. 


Figure  10.  Current  along  x-wall  of  90°  wedge  for  norma 
near  the  shadow  boundary  of  x-wall  by  using 
diffraction  coefficients. 


1 I 
E i 


I 

IV.  CONCLUSION 

The  advantage  of  the  formulation  in  this  report  over  previous 
work  [1]  is  that  the  use  of  the  three  term  series  in  Equations  (19) 
and  (20)  will  permit  one  to  obtain  the  correct  currents  for  all  in- 
cidence angles  whereas  the  previous  formulation  required  special  treat- 
ment of  those  cases  where  the  incident  ray  was  near  grazing  upon  one 
of  the  faces  of  the  wedge.  While  the  three  tern  expansion  in  Equations 
(19)  and  (20)  requires  the  use  of  three  match  points  on  each  face  of 
the  wedge  in  the  GTD  region  instead  of  just  one  (or  two  in  the  case  of 
grazing  incidence)  as  in  the  previous  formulation,  this  is  a small  price 
to  pay  for  a solution  that  works  over  all  incidence  angles. 
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